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WINDOW AVERAGE STATISTICS MODEL FOR POINTING STABILITY 

JITTER ANALYSIS 

BACKGROUND OF THE INVENTION 

5 

[001] The present invention generally relates to pointing stability jitter 
analysis and, more particularly, to a window average statistics model for 
pointing stability jitter analysis. 

[002] Pointing stability jitter analysis is one of the most important analysis 
10 tasks in attitude control subsystem (ACS) design. Pointing stability analysis 
provides the following critical information before the launch of a spacecraft: a 
pointing stability of the spacecraft following a maneuver profile command; a 
pointing robustness at a target in space during a long period of time in a 
mission; a pointing performance of spacecraft attitude jitter with particular ACS 
15 pointing requirements being met or not met. Without an analysis of the pointing 
jitter to make sure the information is within given requirements there would be 
no mission. 

[003] Typically ACS designers will need to run hundreds of simulation 
cases to ensure that all of the simulation cases are within the pointing 

20 requirements after the attitude controller or ACS design is complete. Running 
numerous simulation cases is a time consuming task and places a huge burden 
on the designers who perform the task. A solution is needed that reduces 
substantially and almost eliminates the need for running numerous simulation 
cases by providing a statistics model to almost eliminate the cumbersome 

25 endless simulation runs. The heart of the problem with other people's attempts 
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is the need to run numerous simulation runs. There is a need for a solution that 
cuts down design cycle time, lowers project cost and helps prevent lost 
revenue. 

[004] When analyzing thousands of signals out of hundreds of simulation 
5 runs, a single violation of the pointing stability requirements can cause an ACS 
design, to fail. Since modifications to the ACS are often inconvenient and time 
consuming, failures of this type can be very costly. There is a need to solve the 
pointing violation problem without alterations to the ACS set design. A solution 
is needed to accomplish this by using a statistics model to accurately predict 

1 0 how a particular ACS design will function ahead of time. 

[005] At present, there is no tool for pointing jitter signal analysis that does 
not require altering the ACS design. Without a statistical solution, designers, 
using typical solution methods, would simply overwhelm the ACS computer with 
hundreds, if not thousands, of simulation cases. After collecting the simulation 

15 cases, the designers would hope that signals were all within the necessary 
pointing requirements. 

[006] As can be seen, there is a need for pointing stability jitter analysis in 
an attitude control subsystem that avoids numerous simulation runs. The 
analysis of the information is critical to the launch of a spacecraft. Without a 
20 proper analysis of pointing jitter properties to make sure they are within given 
requirements there can be no mission. There is therefore a need for avoiding 
numerous simulation runs in order to provide the pointing jitter properties that 
are compared with given requirements in determining whether or not a mission 
can proceed. 

25 

SUMMARY OF THE INVENTION 

[007] In one aspect, the present invention provides a method of analyzing 
a signal comprising the steps of: collecting a signal; processing the signal to 
30 produce data; inserting the data from the signal into a statistical model to 
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produce statistics; compiling the statistics using the statistical model to produce 
a probability density function; integrating the probability density function to 
produce a probability distribution function; and determining a pointing accuracy 
from the probability distribution function. 
5 [008] In one aspect of the present invention, a statistics model is disclosed 
for analyzing spacecraft attitude pointing stability in a jitter analysis, the 
statistics model having signals processed to accurately predict the pointing 
stability in flight comprising the steps of: defining an ergodic random process 
statistically and mathematically; creating a window averaging technique to slide 
10 through simulated signals; building the statistics model; loading statistics into a 
statistics metric and creating a probability density function; and integrating the 
PDF to a probability distribution function and reading out a 3-sigma (3-a) 
pointing accuracy against requirements. 

[009] In another aspect of the present invention, a method is disclosed for 
1 5 analyzing spacecraft attitude pointing stability in a jitter analysis by processing a 
limited number of signals to accurately predict the pointing stability in flight, the 
method comprising the steps of: defining an ergodic random process statistically 
and mathematically creating a window averaging technique to slide through the 
simulated signals building the statistics models; placing the statistics into a 
20 histogram and create the Probability Density Function (PDF); and integrating 
the PDF to Probability Distribution function and reading out the 3-a pointing 
accuracy against the requirements. 

[0010] In still another aspect of the present invention, a signal processing 
scheme is disclosed for analyzing spacecraft attitude stability in a jitter analysis. 

25 The scheme processes a limited number of signals to accurately predict a 
pointing stability in flight, the scheme comprising the steps of: collecting a signal 
and processing the signal to produce data. The scheme further comprises 
inserting the data from the signal into a statistical model consisting of statistics 
metric S(t) to produce statistics and deriving the statistics metric using a sliding 

30 window T across the entire signal X(t) one sample at a time. The scheme 



-3- 



PATENT 
091-0195 



further comprises the sliding window T collects S(t) at each window T and S(t) 
is a pointing accuracy. The scheme further comprises compiling the statistics 
using the statistical -model to produce a histogram probability density function 
and integrating the probability density function to produce a probability 
5 distribution function. The scheme further comprises determining a pointing 
accuracy from the probability distribution function where, the pointing accuracy 
is 99.8% of the probability distribution function. 

[0011] These and other features, aspects and advantages of the present 
invention will become better understood with reference to the following 
10 drawings, description and claims. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0012] Figure 1 shows, according to one embodiment of the present 
15 invention, a flow chart of the process of providing the pointing stability jitter 
analysis using a window average statistics model; 

[0013] Figure 2 shows, according to one embodiment of the present 
invention, a random signal, for example, a normal nosy signal with spikes, and a 
window T that slides along the horizontal axis, the horizontal axis is time t, the 
20 window T taking samples of the random signal at X(t), where the vertical axis is 
an amplitude. 

[0014]. Figure 3 shows, according to one embodiment of the present 
invention, a histogram of the random signal of Figure 1 consisting of a 
probability density function after the random signal is put into a statistics metric 
25 S(t). This histogram represents a form of the data derived using the sliding 
window T in Figure 1 . The data is derived by sliding the window T across the 
entire random signal X(t) one sample shift at a time and collecting the data S(T) 
at each window T; 

[0015] Figure 4 shows, according to one embodiment of the present 
30 invention, a probability distribution function curve that results when the data 
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from S(t) as a histogram in Figure 2 is integrated; 

[0016] Figure 5 shows according to one embodiment of the present 
invention, the second method as a statistics metric and a window size of 20 
seconds; and 

5 [0017] Figure 6 shows, according to one embodiment of the present 
invention, the third method as a statistics metric. 

DETAILED DESCRIPTION OF THE INVENTION 

10 [0018] The following detailed description is of the best currently 
contemplated modes of carrying out the invention. The description is not to be 
taken in a limiting sense, but is made merely for the purpose of illustrating the 
general principles of the invention, since the scope of the invention is best 
defined by the appended claims. 

15 [0019] One embodiment solves the point violation problem, in pointing 
stability jitter analysis of the ACS design in spacecraft and satellites. An 
embodiment bypasses the typical analysis concepts by building a statistics 
model around a much smaller set of simulation runs to characterize the overall 
pointing properties. At the end of the analysis, an embodiment will come up a 

20 3-a number for the pointing accuracy to comply with given pointing 
requirements. Using an embodiment, a tremendous amount of project cycle 
time and cost can be saved. Furthermore, an embodiment of the present 
invention provides a more correct analysis than other attempts at solving 
pointing stability analysis or pointing accuracy problems. 

25 [0020] In one embodiment a statistics model is built to analyze a 
spacecraft's attitude pointing stability in a jitter analysis. Processing a limited 
number of signals through an innovative tool set this embodiment will accurately 
predict the pointing stability in flight. The embodiment sets out a signal 
processing scheme which includes: defining an ergodic random process 

30 statistically and mathematically; creating a window averaging technique to slide 
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through simulated signals; building a statistics model, the statistics model is a 
first, second or third method. The next step is placing the statistics into a 
histogram to create the Probability Density Function (PDF1); integrating the 
PDF1 to produce a Probability Distribution Function (PDF2); and reading out a 
5 3-a pointing accuracy then comparing the 3-a pointing accuracy with 
requirements. Sets, for the various methods MATLAB™ software implemented 
routines, are used to process parts of the embodiment in a numerical fashion. 
[0021] A method of analyzing a signal is shown in this embodiment. This 
method comprises the steps of: collecting a signal or a random signal and 

10 processing the signal to produce data. Furthermore, the method includes 
inserting the data from the signal into a statistical model to produce statistics 
and compiling the statistics using the statistical model to produce a probability 
density function. The method of analyzing a signal further includes integrating 
the probability density function to produce a probability distribution function; and 

15 determining a pointing accuracy from the probability distribution function. The 
statistics model uses one of three statistical methods namely the first, second or 
third methods. All three methods use a statistics metric S(t), deriving the 
statistics metric using a sliding window T across the entire signal X(t) one 
sample at a time while the sliding window T collects S(t) at each window T. The 

20 statistics metric S(t) is a pointing accuracy, using the first method, the pointing 
accuracy is a maximum excursion in X(t) from X(to) in window T and S(t) = 
max | X(t+r )-X(t)| where te [0,T]. Using the second method the statistics metric 
S(t), a pointing accuracy, the pointing accuracy is a maximum peak-to-peak in 
window T and S(t) = max (X(t+r )) - min(X(t+r )) where te [0,T]. Using the third 

25 method, the statistics metric S(t),a pointing accuracy, the pointing accuracy is a 
root-mean-square (rms) of X(t) in window T and S(t) = X(t)): X(t+r ) where te 
[0,T]. The probability density function may be a histogram. The pointing 
accuracy may be equal to 99.8% of the probability distribution function. 
[0022] When ACS designers talk about pointing stability in jitter analysis the 
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definition of the signal characteristics and requirements, for example from 
project to project, are often confused. In one embodiment, a correct standard 
for pointing stability jitter analysis is provided. In the past, for example, the 
pointing accuracy was based on a spike 26 as shown in Figure 2. The problem 
5 is that past attempts based the pointing accuracy on spikes that are four times 
as much than the normal signal, normal being the signal without a spike. Past 
attempts would then ignore the spikes without proper bases to do so and more 
importantly would give less accurate results. Therefore, there is a present need 
for a pointing stability jitter analysis that includes the spikes in the signal as part 

10 of the analysis. 

[0023] One embodiment of the present invention is shown in Figure 1 . This 
embodiment of the process is as follows: collect a signal 10, process the signal 
to produce data 12, insert the data from the signal into a statistical model to 
produce statistics 14, compile the statistics using the statistics model to produce 

15 a probability density function 16, integrate the probability density function to 
produce a probability distribution function 18, determine a pointing accuracy 
from the probability distribution function 20. The random signal 22 as shown in 
Figure 2, is real and often seen, for example in some of the rate gyro sensor 
outputs. When the random signal appears over time 24, in seconds, the 

20 pointing stability accuracy could be out of requirements. Instead of attacking 
the problem statistically and mathematically, some designers, in past attempts, 
have spent a lot of time and project money to build software based on control 
laws to correct for the spiky data 26. Such attempts could have been totally 
avoided if the designer knew the solution as shown in the embodiments. This 

25 solution is a statistics model approach. 

[0024] One embodiment first defines a statistics metric in the window T 28 
where T usually comes? from pointing stability requirements. Then, a few 
statistical assumptions are needed to create the statistics model: original signal 
x(t) is a stationary random process; x(t) is ergodic, engodic means statistics 

30 taken in one time history equal statistics taken in different time histories. 
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Hence, sliding the window 30 along x(t) makes statistical sense. A typical 
pointing accuracy requirement states as follows: within time window of 240 sec, 
the spacecraft attitude control error should not exceed 0.02 deg in 3-sigma or 
3a. 

[0025] Once the statistics assumptions are defined, there are three methods 
to specify pointing accuracy S(t) on top of the original signal X(t) in the 
requirement window T: 
First Method: 

Maximum excursion in X(t) from X(to) in window T 

S(t) = m2u\X(t + r)-X(t)\ 
re [0,71 

Second Method: 

Maximum peak-to-peak in window T 

S(t) = max( X(t + r)) - mm(X(t + r)) 
re[0,r] 

[0026] A MATLAB™ code implementing the above equation is simply 
for i=1 : (npts - window + 1 ); 

sample = position( i : (i + window - 1)); 

s(i,1 ) = abs(max(sample) - min(sample)); 

end; 

[0027] Third Method: 

Root-mean-square (rms) of X(t) in window T 

re[0,r] 

A MATLAB™ code implementing the above equation is simply 
for i=1 : (npts - window + 1 ); 

sample = position( i : (i + window - 1 )); 

s(i,1) = std(sample); 

end; 
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[0028] Putting the signal of Figure 1 into statistics metric S(t), we have its 
histogram 32 as in Figure 3. It's derived using a sliding window T across the 
entire signal X(t) one sample shift at a time and collects S(t) at each window T. 
The result S(t) can be put into a histogram 32 in Figure 3. The vertical axis is 
5 probability 34. Figure 3 shows the histogram of the signal in Figure 1 as a 
probability density function 36. Figure 4 shows a curve 40 produced after the 
histogram data in Figure 3 is integrated. Integrating the statistics metric S(t) we 
have a probability distibution function 42. The vertical axis in Figure 4 is 
probability 46 and the horizontal axis is a pointing accuracy 48. 

10 [0029] The histogram of Figure 5, third graph, is integrated 50 to produce 
the Probability Distribution Function 52. This exists where 3-sigma, 6 pointing 
accuracy 54 is, which is also 99.8% of the Probability Distribution Function. 
Figure 5 shows the complete process using the statistics metric of the second 
method and a window size of 20 seconds 56. 

15 [0030] Figure 6 shows the complete process using the third method as the 
statistics metric. Comparing the second method and the third method, we can 
now see that choosing the statistics model is very important for pointing 
accuracy analysis. The end analysis result shows a difference between the 
numbers for pointing accuracy as 1 .2 pointing accuracy 60 in Figure 6, instead 

20 of 6 pointing accuracy 54 in Figure 5. 

[0031] For a spiky signal as shown in Figure 1, we have to use the third 
method to avoid heavy penalty from the spike values. In this embodiment there 
is no need to remove those spikes in the signal with a lot of design and software 
implementation effort because one of the embodiments still meets the pointing 

25 requirements with the spiky signal data included in the analysis. 

[0032] It should be understood, of course, that the foregoing relates to 
preferred embodiments of the invention and that modifications may be made 
without departing from the spirit and scope of the invention as set forth in the 
following claims. 
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